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Abstract 

Analytical solutions to force-free electrodynamics around black holes are funda¬ 
mental for building simple models of accretion disk and jet dynamics. We present a 
(nonexhaustive) classification of complex highest-weight solutions to the force-free 
equations in the near-horizon region of the extremal Kerr black hole. Several classes 
of real magnetically dominated or null solutions, either axisymmetric or nonaxisym- 
metric, are described which admit finite energy and angular momentum with respect 
to the asymptotically flat observer. Subtleties related to the velocity of light surface 
in the near-horizon region are discussed. 
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1 Introduction and outline 

Black holes surrounded by an accreting plasma admit a rich dynamics which is captured 
by the framework of general relativistic magnetohydrodynamics. Under the assumptions 
that the plasma does not feel the Lorentz force and that the matter backreaction is 
negligible, the dynamics reduces to the so-called force-free electrodynamics in the hxed 
background geometry of the black hole [^. This approximation has been numerically 
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shown to be valid close to the poles in certain regimes The force-free eqnations are 

still nonlinear and therefore only few exact analytical solntions are known; see [^|4-11 
for a nearly exhanstive list. In this paper we will continne the recent ongoing effort 
in deriving new analytical solntions which could be used as simple models for accretion 
around rotating black holes. 

Recent observations indicate that near-extremal spinning black holes exist in nature 
12 16 . For such black holes, a throat geometry forms in the vicinity of the horizon. 


This geometry is the so-called near-horizon extremal Kerr region, which admits enhanced 
SL{2,M.) symmetry 17 . Such symmetries are familiar in conformal field theories where 


fields typically fall into highest-weight representations of the conformal group. Building 
highest-weight representations in terms of the electromagnetic gauge field can therefore be 
promoted as a solution generating technique for force-free electrodynamics in the vicinity 
of near-extremal black holes fTol. 


A caveat is that not all the regions of the near-horizon extremal Kerr black hole are 
physical. Indeed, the nonextremal Kerr black hole contains a velocity of light surface, 
dehned as the codimension one surface away from the horizon where the Killing horizon 
generator becomes null, strictly outside of the horizon. The fact that it asymptotes at the 
equator to the horizon in the extremal limit leads to the unphysical presence of a velocity 
of light surface in the near-horizon limit bounded by two polar angles that are symmetric 
with respect to the equator. As a consequence, no global timelike Killing vector exists in 
the near-horizon region and no quantum field vacuum can be defined 18-^ (this was 
also emphasized in [^). So it would not make sense to define a e"*" — e“ plasma in that 
region. This unphysical feature is not present away from the extremal limit. We take 
therefore the point of view that one should build physical solutions in the region between 
the north pole and the northern velocity of light surface only (and similarly in the south 
region). We will call the union of these north and south regions away from the velocity 
of light surfaces the physical near-horizon region. 

Several classes of solutions have been found using the highest-weight technique 
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22,23 . As in each solution generating technique, one distinguishes between the set of 


generated formal solutions and the set of physical solutions. First, highest-weight solu¬ 
tions are generally complex, and only when the current is proportional to its complex 
conjugate up to an arbitrary function could one superpose the real and imaginary parts 
of the solution to get a real solution, as shown in 


22 


Second, many real solutions are 
electrically dominated in the physical near-horizon region. As explained e.g. in 1^, 


electrically dominated force-free solutions are unphysical since there exists a local inertial 
frame for which charged matter moves at a drift velocity higher than the speed of light. 
Null solutions (with electric and magnetic fields of equal magnitude) correspond to the 
borderline case where the charged matter moves at the speed of light. We will encounter 
all such types of solutions but we will seek for magnetically dominated or null solutions. 

Moreover, one should always require that solutions have hnite energy and angular 
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momentum which in the case of highest-weight solutions imposes a bound on their weight. 
Here, there are two possible notions of energy. Either one insists on having hnite energy in 
the near-horizon region, which might be useful for discussing holography in near-horizon 
geometries, or one insists on having a finite and nonvanishing energy with respect to the 
asymptotically flat observer, which we consider as the physical requirement. The latter 
requirement imposes bounds on the weight. For such weights, the variational principle 
should be well-defined in the near-horizon region. 


Our main result is the final list of near-horizon solutions which admit finite and non¬ 
vanishing energy and angular momentum with respect to the asymptotically flat observer 
which we could build from the complex highest-weight solutions. We limited ourselves 
to near-horizon solutions which led to computable energy or angular momentum flux in 
the asymptotically flat region. While we found several classes of axisymmetric or non- 
axisymmetric, magnetically dominated or null solutions, we could not conclude whether 
or not these solutions should be considered physical. Indeed, they are singular at the 
horizon, except one notable null solution found in 22 , and in some cases admit a log¬ 


arithmic divergence at the velocity of light surface. One can argue that the velocity of 
light surface of the near-horizon region should have order one corrections upon gluing 
the asymptotically flat region. If this is the case, this logarithmic divergence might be 
unphysical. This issue could only be totally settled by considering the extension of the 
near-horizon solutions in the asymptotically flat geometry, which is a difficult analytical 
problem. 

Our paper is organized as follows. After a brief review of the near-horizon extremal 
Kerr geometry, we first present an extension of the formalism of Euler potentials 


25 to 


describe nonaxisymmetric solutions in a canonical form in Sec. which will be instru¬ 
mental in expressing the solutions. We then discuss the S'L(2,]R) invariant solutions in 
detail in order to get some intuition. Our main analysis takes place in Sec. First, we 
define a S'L(2,M) covariant basis of 1- and 2-forms to express all the physical quantities 
in a covariant fashion. We then construct the highest-weight representations of the vector 
potential, held strength and current. Second, in Sec. |5.3|, we analyse the energy and 


angular momentum huxes with respect to both an observer in the near-horizon region 
and an asymptotically hat observer. We are interested in solutions with hnite and non¬ 
vanishing asymptotically hat energy and angular momentum that are determined from 


the near-horizon region. This analysis, according to 26 , provides bounds on the weights 


of the solutions. Third, we show that in order to obtain highest-weight representations 
obeying the force-free equations we only need to solve three coupled nonlinear ordinary 
diherential equations. We classify the types of solutions using several criteria (electro¬ 
magnetic type, admitting or not admitting descendants, etc) and partially reduce the 
equations to a list of independent solutions. We however identify one class of two coupled 
nonlinear equations that we were not able to reduce, so our classihcation, even though it 
extends the existing literature, is not complete.We then list all the independent complex 
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highest-weight solutions which we obtained. In Sec. we hrst build real solutions, from 
the list of complex solutions, which are magnetically dominated or null. This preliminary 
result leads to a list of near-horizon solutions in Sec. 6.2 By further imposing the bounds 


derived in Sec. 5^ to have hnite and nonvanishing energy and angular momentum which 
respect to the asymptotically flat frame, we get the hnal list of potentially physical near¬ 
horizon solutions in Sec. H The highest-weight solutions are expressed in terms of hve 
distinct functions of the polar coordinate. Out of these hve functions, three solve lin¬ 
ear ordinary differential equations (ODEs) and two solve nonlinear ODEs with boundary 
conditions hxed by regularity requirements as we discuss in Appendix We numerically 
solved explicitly all three linear ODEs. For completeness and future reference, the details 
of all highest-weight solutions found are listed in Appendix 

Note added in v3 After the publication of (2^, we noted two algebraic mistakes in 


Eq. (97) and (101) in v2 which we corrected in this version, in Eq. (97) and (102). These 
corrections modify the admissible weights and the list of potentially physical solutions. 
Moreover, complete bounds on the weights require the analysis of angular momentum as 
well, as done in 26 . Sec. 5.3 and 6.3 have been rewritten and the introduction and 
conclusion have been adapted accordingly. We also commented further upon regularity 
at the future Poincare horizon. 


2 Near-horizon extremal Kerr geometry 


The near-horizon limit of the near-extremal Kerr black hole has been explicitly written in 
many references and we refer the reader to those 


17 27 2') 


It is important to emphasize 
that depending on how the near-extremality parameter J — is scaled upon taking 
the near-horizon limit, two distinct coordinate systems might result in the near-horizon 
extremal Kerr geometry, namely Poincare and black hole coordinates. Such coordinate 
systems do not cover the entire near-horizon spacetime and can be extended into global 
coordinates. It is therefore of interest to write a formalism which does not depend upon 
the choice of coordinate system and which can be easily specialized for each case. 

The four-dimensional near-horizon extremal Kerr metric can be written in a way that 
makes its S'L(2,M) x U{1) symmetries manifest but without choosing a particular coor¬ 
dinate system as 


ds^ = p(0) [dslas, + + 'i\d) [dT + kef] . (1) 

Here T ~ T -|- 27r is the azimuthal angle, 6 G [0, tt] the polar angle, ds\^g^ is the unit 
metric on two-dimensional anti-de Sitter spacetime AdS 2 and 0 is a left-invariant one- 
form on AdS 2 with norm —1. The geometry depends upon the two functions 7 , P of the 
polar angle and upon the constant k. We will keep k explicitly here in order to more 
easily allow for future generalizations to other near-horizon geometries such as e.g. the 
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one coming from the Kerr-Newman black hole which has the form Q with k ^ 1. For 
the extremal Kerr near-horizon geometry, we have 


k = 1, 


r(0) 


1 -I- cos^ 6 


7 ( 0 ) 


2 sin0 
1 -I- cos^ 6 


( 2 ) 


where J is the angnlar momentnm. The velocity of light snrface is located in the range 
[0*, TT — 0*] where 0* = arcsin(\/3 — 1). We will mostly consider the physical near-horizon 
region defined as [0, 0*] U [vr — 0*, tt]. 

We denote the symmetry generators of S'L(2,]R) as Hq,H± and the U{1) symmetry 
generator as Qo = dq,. For fntnre nse, we also define the rescaled generator Qo = q^.q^ ■ 

When considered as a form after using the metric, it has the convenient expression Qo = 
dT -|- kO (we denote forms and vectors by the same symbol and they are distinguished by 
the context). 

The generators obey the S'L(2,M) x U{1) commutation relations 


[Ho, H^] = [H+, H.] = 2Ho, [Qo, H^] = 0, [Qo, Ho] = 0. (3) 


A physical feature of this geometry is the absence of globally timelike Killing vector. 
Here, we complement the comments given in the Introduction. For a nonextremal black 
hole in comoving coordinates the Killing horizon generator is null at the horizon and 
timelike just outside the horizon and it becomes null away from the horizon at the velocity 
of light surface beyond which it becomes spacelike. At extremality and at the equator 
the velocity of light surface asymptotes towards the horizon and therefore upon taking 
the near-horizon limit, the Killing generator might not remain timelike everywhere. The 
region where the Killing generator remains timelike represents a region where the physics 
is clearly related to the asymptotically flat region. On the contrary, the region where 
it becomes spacelike is rather a special feature of the near-horizon extremal limit which 
disappears in the asymptotically flat region. 

For further use, we define the highest-weight scalar ^{h,q) of weight h and charge q as 


H+^{h,q) — 0, 

Ho^(h,q) — h<h(h,g), 

Qo‘^ih,q) = 

dg^ = 0. 

(4) 

We denote 

III 

III 

ye 


(5) 

Since the definition is 

linear it readily follows that 




^{h,q) - 



(6) 


We will now make contact in the upcoming sections with the three special coordinate 
systems: Poincare coordinates, global coordinates and black hole coordinates. 
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2.1 Poincare coordinates 

We denote Poincare coordinates as {t,r,6,(j)). We have \1/ = 0. The metric is 


We dehne 


ds^ = r 


dr^ 

+ d6^ — r'^dt^ + 7 ^ [dc/* + krdt\ 


H+ 

Ho 

H_ 

Qq 


V2dt, 
tdt - rdr, 



(t^ + \)dt - trdr - -d, 


We have Qo = d0 + krdt. The highest-weight scalars are given by 


$ = -, 
r 


A = 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


( 12 ) 


2.2 Global coordinates 


We denote global coordinates as {T,y,9,i^). We have '1 = 9 ?. The metric is 

^ dy^ 


ds^ = T 


1 + 1/2 


+ d9^ — (1 + yjdr + Y [d^ + kydr]" 


We define 


77+ 

^IT 

-7 ( 

'-ydr + i{l + y‘^)dy - kd^, 


77o 

= idr 


H_ 

p-ir 

-i ( 

'-ydr - i(l + y‘^)dy - kd^, 


Qo 

= d^. 



We have Qo = dip + kydr. The highest-weight scalars are 

_ p- 2 r 

$ = \ = g++fearctany_ 

+T7 


(13) 


(14) 

(15) 

(16) 
(17) 


( 18 ) 
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2.3 Black hole coordinates 

We denote black hole coordinates as We have The metric is 


ds^ = r 


dY^ 


-1 + ^2 


+ de^ - (-1 + Y^)dT^ + Y [d'lp + kYdT] 


( 19 ) 


We dehne 


n-T 


H+ = 

J-^iYdr + (-1 + Y^)dY - kd^), 
v-1 

(20) 

Ho = 

dx 

(21) 

H_ = 

/ 1 1 ^ l + Y’^)dY kdY, 

V—1 + y 

(22) 

Qo = 

d^. 

(23) 

+ kYdT. 

The highest-weight scalars are 



^ = V2 


x/-l + d"2’ 


^ z/c arctanh y 


(24) 


3 Canonical Enler potentials 


The use of differential geometry in the covariant formulation of force-free electrodynamics 
was elegantly motivated and developed in 30 . We will follow their conventions. Maxwell’s 
equations are 

dF = 0, d-kF = -kJ (25) 


where F is the electromagnetic held 2-form, J is the current 1-form and * is the Hodge 
star operator. The force-free condition is expressed as 


JA4^F = 0. (26) 

Force-free electromagnetic helds are degenerate, i.e., det(F) = 0. This implies that the 
electromagnetic held can be written as 


F = d(j)i A (i02 


(27) 


where and 02 are the so-called Euler potentials which are determined up to a held 
redehnition of (0i,02) of unit Jacobian. 

We only consider spacetimes such that the metric admits a block diagonal form into 
a so-called two-dimensional Lorentzian toroidal part and a two-dimensional Euclidean 
poloidal part. The Kerr metric admits this decomposition and it follows that the near¬ 
horizon extremal metric also does. 












Let us now concentrate on the near-horizon extremal Kerr in Poincare coordinates. We 
anticipate that highest-weight solutions are constrained by iL+K = 0 which is equivalent to 
requiring stationarity of the held strength. Canonical Euler potentials have been derived 
in this context for stationary and axisymmetric conhgurations 25 . Here, as a primer 


for describing the physical properties of highest-weight Poincare solutions we hnd useful 
to hrst derive (complex) canonical Euler potentials for stationary but nonaxisymmetric 
conhgurations which are (complex) eigenstates of 5^, QqF = iqF. 

The metric in Poincare coordinates can be decomposed into the toroidal part spanned 
by (t, 0) with volume form er = Pyrdf A dcf), and the poloidal part spanned by (r, 6) with 


volume form = -dr A d9. 

r 


We have e = A 


-kC — —6 , 


'ke — 6 


3.1 Stationary and axisymmetric case 


Let us hrst summarize the stationary and axisymmetric case as analyzed in 25 and 
reviewed in [^. There is no toroidal electric held (and therefore no components of 
the held strength proportional to dt A dcf)) for axisymmetric conhgurations as a simple 
consequence of Faraday’s law. We distinguish three scenarios: 

generic case ig^F ^ 0. One can choose 

01 = ?/>(r, 6 '), 02 = 0 + 02 (a^) - ^(^)L (28) 


The polar current I (r, 9) is dehned as 

* {dip A dip2) = ^ ^ e^. (29) 

It is equal to the electric current with respect to time t howing in the upward 
direction through the loop of revolution dehned by the poloidal point {r,9). Note 
that this interpretation breaks down beyond the velocity of light surface where dt 
is spacelike. The force-free equations imply that I = I{ip{r,9)). We therefore have 

Ht a df) 

F = dip A {d(p — VL{ip)dt)+I{ip) - - —. (30) 

In particular, if VL{ip) = 0, there is no electric held, Iq^F = 0. 


No poloidal magnetic field ia^F = 0, ig^F ^ 0. One can choose instead 

0i = x(a^), <p2 = tFX2{r,9). (31) 

We then dehne the polar current as '^{dx A dxi) = dPpPLpp which has the same 

interpretation as above. The force-free equations imply / = I{x{r,9)). The corre¬ 
sponding held strength takes the form 

^ ^dr A d9 , , 

F = dxAdt + I{x) - (32) 
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Only toroidal magnetic field ig^F = 0, ig^F = 0. In that case, 

fJv A df) 

0i=x(r,0), 02 = X2(r,0), F = I{x) -(33) 

There is no electric field and no poloidal magnetic field. 

3.2 Stationary and d^-eigenvalue case 

Let us now consider a complex force-free field strength which is stationary, Cg^F = 0, and 
which is an iq d<ji-eigenvalue, Cg^F = iqF. Stationarity implies 

0 = Cg^F = dig^F = d (-5*02^01 -I- 0t0id02) (34) 

where we used Cartan’s formula, Bianchi’s identity and the degeneracy of F. By Poincare’s 
lemma, there exists a function / = /(0i, 02 ) such that: 

- (9*02^01 + 0t0i(i02 = df. (35) 

We distinguish here two cases (i) ig^F = df = 0 which implies that both Euler potentials 
are time independent and (ii) ig^F = df ^ 0 to which we now turn our attention. Euler 
potentials are defined up to the following arbitrariness: we may choose any other pair 
of potentials (0i,02), leaving the electromagnetic 2-form invariant, provided the map 
(01,02) —t (01,02) has unit Jacobian determinant. Using this freedom, we choose 0i = 
—/. Let us check the existence of 02(0i,02)- The Jacobian of the transformation reads 
as 


1 = 


001 002 002 001 


df 002 ^ 0 / 002 


001 002 002 001 


(36) 


which is a first order partial differential equation (PDE) for 02 ( 0 i, 02 ) and can be inte¬ 
grated with respect to 02 if ^ 7 ^ 0 or with respect to 0i if ^ 7 ^ 0. With this new pair 
of Euler potentials, Eq. (35) becomes 

- 0102^01 + 0t01#2 = -d01 


from which we read off the conditions 

0i01 = 0, 0102 = 1 

whose solutions are 

A = Xi(u 6 ', 0 ), 02 = t + X 2 (r, 6 ', 0 ). 


(37) 


(38) 


(39) 
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Finally, merging cases (i) and (ii) and dropping tildes, Euler potentials for stationary 
solutions can be fixed to 


(f>i = (j )2 = et +X2{r,9,(j)) (40) 

where e = 1 if ig^F ^ 0 and e = 0 if ig^F = 0. 

Let us now turn our attention to the second condition Cg^F = iqF. We have 

0 = dig^F - iqF 
= dig^{d(f)i A ^02) — iqd(pi A d(j )2 

= d [(i6»^d0i)d02 - - iqd(t)i A d02 (41) 

= d [{d^4>i - iq(t>i)d(j)2 - {d^4>2)d4>i] 

= d [{d4>xi - 'iqXi){.^dt + dx2) - {,d^X2)dxi] 

where we used Bianchi identity in the hrst step and stationarity in the last one. Let us 
hrst discuss the case e = 1. Since Xi,X 2 have no time dependence, from the identity 

0 = d{d^xi - iqXi) A {dt + dx2) - d{d^X2) A dxi (42) 

we infer that 


< 90 X 1 - iqxi = const, 8 ^X 2 = «;(Xi) (43) 

where k{xi) is an arbitrary function of the Euler potential xi ^cid where the arbitrary 
constant can be set to zero by shifting xi- 
From the hrst differential equation we have 

Xi(r,0,0) = e*‘?‘^Xi(r,0). (44) 

From the second differential equation, we infer 


X2{r,0,il>)= I + X2(r,9) 


(45) 


In conclusion, dropping the tildes for simplicity, the Euler potentials in the case ig^F 7 ^ 0 
can be taken as 


01 = e*'''^Xi(r,0), 02 = t + X2(r,0) + / K(e*'?'^'xi(r, 0))ci0' 


(46) 


Let us compute the held strength. We dehne h{r, 9, 0) = K{e'^i^'xi{r,9))d(j)', then 

11 f n d)h ^ dh dh ^ ^ 

dh{r,9,<i)) =—dr+—d9+—d<i) (47) 
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where 


dh 

dr 




d(j) 


dh 

T4, = 


The exterior derivative of 02 is 

d(j)2 = dt + dx2 + K(0i)d0 + 
and the field strength takes the following form 


5«)(01 

8 h 


1 dxi 


( 49 ) 


(50) 


F = d(j)i A {dt + dx2) + e 




- iQXi j 


dXi A dcj). (51) 


Let ns now retnrn to the case e = 0 {idtF = 0). We restart from (41). By Poincare’s 
lemma there exists a fnnction / snch that 


[{d^Xi - iqXi){dX2) - {d4>X2)dxi] = df. 

li df = 0 we find directly 

d4>X2 = 0 , d^xi - iQXi = 0 

and we find the Euler potentials in the case dtF = 0 with df = 0, 


<t>i= e“'*Xi(r,0), h = X 2 (r,ff). 


(52) 

(53) 

(54) 

(55) 


In that case, the held strength is 

F = (^dxi A dx2 - iqXidX2 A df'^. 

If df 7 ^ 0, one can use again the ambiguity in the dehnition of Euler potentials to choose 
Xi = -/• Then 

d4>X2 = 1, d^xi - iqXi = 0 (56) 

and we hnd the Euler potentials 0i = e*'^'^Xi(r, 6*), 02 = 4> + X 2 {'f', 0). Since these potentials 
generalize (28) when hi = 0, we hnd it convenient to align the notations so that hnally we 

(57) 


get in the case dtF = 0 with df ^ 0, 


0^ = (r, 0), 02 = 0 + 02 {r, 0 ). 


In that case, the held strength is 

F = (^dfj A d02 + {df) — iqfjdfj2) A df ^. 


(58) 
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4 Maximally symmetric solutions 

Before considering highest-weight solntions, let us hrst obtain all SL{2, M) x f/(l) solutions 
to force-free electrodynamics. This analysis completes the one of [lo] . 

In the gauge Aq = 0, symmetries imply that the gauge potential has the form 

A = Aq{6)Qo. (59) 


This form is expressed only in terms of the metric and the U{1) Killing vector and is 
therefore S'L(2,M) x f/(l) invariant. Reality of the gauge potential requires that ^ 0 ( 6 ^) 
be real. The current is given by 


J 


1 

T 


a (^) + 


Qo- 


(60) 


The force-free condition is 


0 = J A 


deAo 

T 





+ k'^lAo 


dt A dr A dO 


(61) 


Therefore the only force-free solution with nontrivial current is Aq = —EQ/k constant 
which leads to 


A = -^Qq, F = -^dQQ, J = Eq^Qq 
k k f 

The solution obeys 


(62) 


*{F A *F) = -jU-- 


pfiiA _ AAl 
Y2 


(63) 


In the physical region 1 — > 0 where the Killing vector = \/2dt is timelike, the 

potential is electric with respect to the Killing time t. It is therefore electrically dominated 
and the plasma which underlies the current therefore moves at ultrarelativistic speed; 
see e.g. 24,31 . We will therefore not consider this solution further. Note that these 
conclusions did not depend upon the Kerr functions ([^. It is a consequence of symmetry 
alone. 


The other solution to (61) has no current J and therefore obeys free Maxwell’s equa¬ 
tions. It is given by 


^ — Aq{6)Qq — Mqcos[9q — k j d9'''y{6')\QQ^ 
= Mq cos[ 6^0 -l- 2 arctan cos 9]Qq, 


( 64 ) 
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where 9q is a phase and Mq is a constant magnitnde. In the last line we specialized 
to the Kerr case ([^. The held strength is regnlar at the north and sonth poles since 
^o(O) ~ ^o(^) ~ 0 - The electromagnetic invariant is given by 

* (F A *F) = (-(deA„f + (65) 

We reqnire that the held strength be magnetically dominated [^{F A -kF) < 0] in the 
physical region, ontside the velocity of light snrface. It tnrns ont that it is possible to do 
so npon choosing the phase in the range 

— 2 arctan[Y^2 a/3 — 3] + — < 6^0 < 2 arctan['\/2 a/3 — 3] — — ( 66 ) 

as clear from Fig. 


*(F A *F)(0) 



00 = -0.41 
00 = 0 
00 = +0.41 


Figure 1: ^{F A-kF) is negative outside the veloeity of light surface, whose boundaries 
are highlighted, if the phase Oq is in the range ( 66 ). 


It is interesting that a maximally symmetric magnetic-type solntion exists in the near¬ 
horizon region. However, it is not snstained by matter helds so it is nnclear how such 
a vacuum solution would be compatible after backreaction with “no hair theorems” for 
stationary axisymmetric black hole solutions to Einstein-Maxwell theory. 


5 Highest-weight solutions 


Our objective is to classify the (complex) highest-weight solutions of force-free electro¬ 
dynamics in the near-horizon region of extremal Kerr, thereby completing the analysis 
of 22 . The hrst step in this program is to choose a basis for SL{2, M) x f/(l). There is a 
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large automorphism group of this algebra which allows us to change basis while preserving 
the commutation relations ([^. Continuous automorphisms are parametrized by complex 


(a, /3, 7 , 5) and given by 

H+ ^ e-^ [(1 + a/3YH+ + 2/3(1 + a/3)Ho + P^H_] , (67) 

Hq —)■ a{l + af3)H^ + {1 + 2af3)Ho + f3H_, ( 68 ) 

H_ [a^H+ + 2aHo + H_] , (69) 

Qo SQo- (70) 


The rescaling of Qo parametrized by 6 can be absorbed into a rescaling of the angle </> and 
we do not consider it any further. Note that when a = e'*' = /3^, 5 = 1 and in the 

limit /3 —)■ 0 one hnds the discrete automorphism given by 


H+ H 




Hq —)■ —Hq, Qo ^ Q 


0 - 


(71) 


Given a S'L(2,M) x f/(l) basis, we look for solutions to the force-free equations (25)-(26) 
satisfying the highest-weight conditions 


Ch+F — 0 
FhqF = hF 
FqoF = iqF 


(72) 


where h is real or complex and q is an integer. The choice of automorphism labeled by 7 is 


irrelevant: the resulting solutions to (72) will be identical. Also, it has been shown that the 
solutions for a given set H±, Hq and for a set transformed with real a,/3 can be mapped 
into each other with a real change of coordinates, which is an isometry of the metric. 
Therefore, one might just choose a particular basis and ignore such maps. However, 
the automorphisms labeled by complex a, f3 map a given solution to another nontrivial 
solution since they cannot be mapped into each other via a change of coordinates, as 
noted in 22 . In summary, the automorphism group S'L(2,C) = SL(2,M) x S'L(2,M) 
gets quotiented out by the real transformations S'L(2,M) and the remaining real scaling 
transformations 7 G 77(1) and the resulting nontrivial choice of SL(2,'K) basis is labeled 
by an element (q:,/ 3) of S'L(2, M)/77(l). 

Since the force-free equations are nonlinear, several classes of solutions will appear and 
for each class, one could consider an extension of these solutions using the SL{2, ]R)/t7(l) 
complex automorphisms. Let us now £x a SL{2, M) x 77(1) basis and classify the solutions. 


5.1 SL{2,M.) covariant basis of forms 

Since we would like to keep S'L(2,M) covariance manifest, we need a SL{2,'R) covariant 
dehnition of a basis of spacetime 1 - and 2 -forms that we will use to express all physical 
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quantities. Also, since we describe highest-weight representations, it is computationally 
advantageous even though it is not necessary to choose a highest-weight basis. 

A basis of 1-forms p®, i = 1,... 4 of highest weight and weight 0, = 0, 

is the following 


/i 

/i 


1 

2 

3 

4 


de 

wy 

(i<* 


* ’ 


(73) 

(74) 

(75) 

(76) 


where <h is the highest-weight scalar of weight 1 dehned in (|^. In Poincare, global and 
black hole coordinates we have 


= 


/i = 


= 


/i = 


, , idy dY 

rdt = ydr + ^—77 = Y dT + 


1 + 2/2 - ■ l-y2’ 

d(j) = dip - =di} - k ^ , 


dO 

wr 


1+2/2 

ydy 


dr 

— = 1 , 2 

r 1 + 1 / 


= dT + 


ydy 

i-y2' 


(77) 

(78) 

(79) 

(80) 


A basis of 2-forms w\ i = 1,... 6 of highest weight and weight 1, Ch+w^ = ChqW'^ — w'^ = 0, 
is given by 


A 

(81) 

nP' = A yP, 

(82) 

nP = $/c7^(5o a 

(83) 

= $dQo, 

(84) 

= $Qo A yP, 

(85) 

$2 

A Qo- 

(86) 


This choice of basis is motivated by the Hodge duality properties 

-kW^ = —W^, -kVO^ = —W^, kW^ = —W‘^, (87) 
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and by the properties under the action of the total derivative 


dw^ = 0, Vz 7 ^ 3, 


^3 rj 

-kdw = —pz—Hz 

V2r^ 


( 88 ) 


We could trade the basis elements for another pair with -kw"^ = —w'^ and dw^ = 

0 = dw^ but it would involve double integrals of 7 so we prefer not to use such a basis. 


The normalization of <h in (18)-(24) was chosen so that the duality transformations (87) 
are obeyed. 

Since the basis of 1- and 2-forms is real in Poincare and black hole coordinates, these 
bases are well suited for studying real solutions as we will describe later on in Sec. 


5.2 Vector potential, field strength and current 

We expand the electromagnetic field tensor in the basis w* 


where we chose 


C-H+W'' = 0 , = w\ 


After some simple algebra, the system (72) becomes 


r H+Fi = 0, 

I HoF, = {h- 1 )F„ 

[QoFi = iqFp 


(89) 


(90) 


(91) 


where H+Fi = H^d^Fi, etc. 

We therefore reduced the problem to six decoupled scalar equations of the form (|^ 
whose solutions were provided in ([^. The most general expression for the electromagnetic 
held F is therefore given by 


F(h,q) — ^{h-l,q)fi{d)w'‘ (92) 

where is the highest-weight scalar of weight h — 1 and charge q. The Bianchi 

identities dF = 0 reduce the number of independent functions fi{0) from six to three. 
Then we must hx the remaining three functions by using the nonlinear equations of 
motion. 

The most general vector potential g) generating F(^h,q} will also contain three func¬ 
tions of 6. We consider the following vector potential 

^{h,q) = *h(/i,g)ai(6')/z*. (93) 
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Later on we will work in the gange 04 ( 0 ) = 0 bnt we keep this fnnction arbitrary in this 
section. The held strength fnnctions fi{0) are then determined by 


' fi{d) = ha^ie) 

f 2 id) = -ikqasiO) + 7 i-a[{9) + ka'^iO)), 

^ 7 ^ 3 ( 6 *) = iqasiO) - 

kfiiO) = (1 — h)ai{6) + k {ha 2 { 0 ) — iga^iO )), 
= -ha2{0) + iqa4{9), 

JeiS) =iqai{9). 


Of course, the Bianchi identities are then satished. Also, given an arbitrary set of held 
strength functions which obey the Bianchi identity, one can invert the system in a given 
gauge and solve for the gauge potential functions. Any highest-weight held strength can 
then be expressed from a highest-weight gauge potential. 

From the held strength, one can obtain the current as 


J{h,q) ^{h,q)ji{9)fi 

where 


(95) 


= 7 /' + iqU {k^i^ - 1 ) + ^7^ N/s + (^ - /4 + (h - l)/ 5 ] , 

rj 2 = {h-l)f5 + ikqf6-k'j{k'jf4-fs) , 

rj 3 =-{h-l)fi-ikq{f 2 +fs), 

.7^rj4 = 7/{ - (^V /4 +/s) • 


5.3 Energy and angular momentum flux and variational princi¬ 
ple 

Let us pause to compute the energy and angular momentum hux. In Poincare coordinates 
we dehne the energy with respect to dt, the angular momentum with respect to —d^ and 
the boundary is at r —)■ cx) and the horizon is at r = 0. The hux density of energy and 
angular momentum per unit time per solid angle is given by 

S = ^T^^,{dtyn^ = E{9)r^-^^e^^^'^ (97) 

j = -^Ty{dy^n^ = (98) 

where is the unit normal, 7 ^ 1 . is the induced metric on constant r surfaces and 
is the electromagnetic stress tensor. For simplicity we only discuss complex helds. Real 
helds will be discussed in Section 6.1 The functions E{9), J{9) are explicitly given by 

(99) 

7 / 

■ qa3(,9)\ 


E{9) = hay9)a[{9) + iqay9) ( {h - l)k-fay9) + h- — ^—a2{9) 


J{9) = ha2,{9) ^— 02 ( 6 ^) -|- — iqai{9)'j(^{h — l)ai( 6 ') — hka2{9) 


( 100 ) 
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Requiring no energy and angular momentum flux at the boundary of the near-horizon 
region is, assuming E{6) 7 ^ 0, equivalent to requiring Reh >10 The energy flux however 
then diverges at the Poincare horizon. 

Instead, let us define the energy with respect to the physical time from the point of 
view of the asymptotically flat observer. Let us remember how the near-horizon limit is 
taken. The outside asymptotically flat time, radius and polar angle scale as tout = t/X, 
Tout = M + Xr, (pout = 0 + Xloxtt/X where A —)■ 0. Here M is the extreme Kerr mass and 
VLext = ^ is the extremal angular velocity. Then + VLextd'^^ = Xdt. Since the physical 
energy 8out and angular momentum Jout are associated with 5°“* and the physical 

deviation of energy flux density with respect to comoving flux density per asymptotically 
flat time unit and solid angle can be expressed as 


Jo.. = y=t^r'‘j“‘>(8r‘ + s2«(ar)''d°“‘) 


( 101 ) 


where the right-hand side is a function of {tout, 'Tout, (pout)- We require that the stress-tensor 
admits a well-defined near-horizon limit. In order to reach stationarity in the near-horizon 
limit, we also impose that the stress-tensor only depends upon {rout, (pout — X^exttout) in the 
near-horizon limit, i.e. the held co-rotates with the black hole at the horizon. Now when 
substituting Tout as a function of r, all r factors come equipped with a A factor as well. 


From the near-horizon scaling in Eq. (97) we can deduce the scaling 

C _ O 'T r-j )i2-2/ic 
^out ^ ^extOout — ^ 


( 102 ) 


Following the same reasoning, the angular momentum flux density scales as 




(103) 


The scaling limits (102)-(103) agree with 


Physical solutions should have a finite (noninflnite) physical energy and angular mo¬ 
mentum flux in the limit A —)■ 0. This requires Re(h) < 1 or E{9) = 0, \/9. This condition 
also follows from requiring flniteness of the energy flux density ( [9^ at the Poincare horizon 
in the near-horizon region. Angular momentum should also be finite in the asymptotically 


flat region. Comparing Eq. (103) with Eq. (98), this is equivalent to requiring flniteness 
at the Poincare horizon in the near-horizon region. This imposes that Re(h) < | or 
J{9) = 0, V0. 


Let us discuss the different cases similarly to 26 


(a) Re(h) < one has Sout = Jout = 0 and there is no energy flux extraction from the 
near-horizon region; 


(b) Re{h) = ^, one finds jut = ^extJout with jout - J] 


^This bound can also be obtained in global and black hole coordinates since these coordinate systems 
admit the same falloff at the spatial boundary. 
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(c) ^ < Re{h) < 1 and J{6) = 0, one has Sout = ^extJout with (and hence £out) 
nndetermined; 

(d) Re(h) = 1 and J{9) = 0 one has Sout — S^extJout — S with Jout nndetermined; 

(e) Re(h) > 1 and J{6) = E{6) = Q and Sout and Jout are not determined from the 
near-horizon region information. 


Here we allow for complex h, which indicates an oscillatory behavior at the spatial bonnd- 
ary. In the case (b) and (c) the energy flnx satnrates the snperradiant bonnd. This indi¬ 
cates that the field exactly corotates with the black hole. It can only happen for photons 
orbiting at the Killing horizon which classically never escape it. This is a featnre of the 
ideal extremal solntion. In the case (d) one has < 0. Indeed, as clearly explained in 30 , 
npon changing the mass and angnlar momentnm of the black hole as 5M = —SSout and 
6J = —5 Jout with ontgoing radiation we mnst obey 5M — Qh6J = -^5A > 0. Therefore 
Sout < S^extJout- This physically corresponds to the fact that an ontflow of energy nec¬ 
essarily comes accompanied with an ontflow of angnlar momentnm. This ineqnality can 
also be obtained from the nnll energy condition ^ 0 nsing the Killing horizon 

generator x = 9t + Sind,/,. 

We will be interested in solntions which admit an energy flnx in the asymptotic region 
which implies Re(h) > We will limit onrselves to near-horizon solntions which bring 
information abont the asymptotically flat region, which implies Re(h) < 1. Onr definition 


of potentially physical solntions in Section 6.3 will only refer to solntions with 


- < Re(h) < 1 and 


which covers the cases (b), (c) and (d). 


Re{h) - - 


J{0) = 0, 


(104) 


The bonnd (104) also has a special meaning in the variational principle. Indeed, the 
action of force-free electrodynamics is given by 

= — - [ dj A d(j)2 A -k^dj A d(j)2) (105) 


whose variation gives a bonndary term Q'^{d^x)fj, with 

0 ^ = y/^F^^^{5(t)2dy(t)i - 


(106) 


where is the held strength. It is easy to check that for a 

highest-weight solntion one has 

0" ~ (107) 


Therefore the case Re h = 1 is precisely the one for which there is symplectic hnx, and 
additional care is reqnired to dehne the variational principle. The variational principle is 
otherwise well-dehned in the range (104). 
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5.4 Solving the force-free condition 


For the highest-weight ansatz, the force-free condition (26) is equivalent to the following 
relations between /j’s and jVs: 


k{f2 + /3)j3 + (/s + A;^7^/4)j4 — /eji, 

/2j3 + fc7V4j4 = hj2, 

-Wfijl + (/5 + /cV/4)j2 = fljs, 

J2jl -k{f2 + f3)j2 = fljA- 

Now, we see that we can recast this system of equations in the following form 


(108) 


A 

-/el 

X 


o' 

-/il 

B 

y. 


0 


where 


and 


A = 


k{f 2 + fs) fb + ^^7^/4 
/2 


= aA^a, 


a = 


X = 


0 1 

-1 0 


y = 


(109) 


( 110 ) 


( 111 ) 


The linear system (|109|) has nontrivial solutions if its determinant is vanishing, i.e., 

4 -/el 


0 = det 


-/il 


— [det(A) -|- /i/e] , 


( 112 ) 


which turns out to be equivalent to the degeneracy condition of the held strength F. For 
dehniteness, we work in the gauge 


04 = 0 Vh and 03 = 04 = 0 for h = 0. 
The degeneracy condition can then be written as 


{h — l)aia2 — ha)02 -|- iq 


0403 

7 


= 0 . 


(113) 


(114) 


Substituting (94) and (96) in (108) one obtains in total three nonlinear coupled ODEs 
in terms of the gauge potential functions 04 , 02 , 03 . The hrst equation is (114). The two 
other equations are lengthy and unenlightening. Since an interested reader can easily 
reproduce them we will omit them here. 

Given the difficulty of these equations, it is useful to organize the solutions by hrst 
classifying the solutions in terms of held strength functions only. After analysis, we could 
branch these nonlinear equations into seven independent and complete (still nonlinear) 
subcases where they can be solved: 
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1- /i — /2 — /s — /4 — /e — 0, /s 7^ 0, j2 — Ja — 0, 

2. /i = /2 = /s = /e = 0, U 7^ 0, Ja = 0, ji = kj2 + 


fc/ 47 ^ 


3. fl = f2=fA = h = 0, /s 7^ 0, J2 = 0, J3 = -ft, 


4. fi = 0, /2 ^ 0, /s = 


_ fc^7^/3/4 _j_ kf3j2 — /6i2-fc7^/4i4 


/2 


/2 


J3 = 


/2 


5. /i = A = A = 0, A A 0, Ji = ^ 

6 . A = A = A = 0, A A 0, J2 = 0, ji = 


7 . A 7 ^ 0 , A = J 3 = 


_ /5i2+fc7^/4(-il+fci2) ■ _ -kf3j2+f2ijl-kj2) 


fl 


3 a = 


h 


Substituting (94) and (96) one obtains the nonlinear equations in terms of the gauge 
potential functions a*, i = 1, 2, 3. We were able to fully solve the cases 1 to 6. However, 
the general solution in case 7 eluded us. The resulting two ODEs are strongly nonlinear 


and it is not clear whether we obtained all possible solutions. We will present in Sec. |5.7 
all solutions found. 


5.5 Electromagnetic types 

In order to better present the solutions, it is useful to hrst present some criteria for distin¬ 
guishing magnetically dominated solutions from null solutions and electrically dominated 
ones. 

We dehne the (not normalized) electric and magnetic helds with respect to the vector 
as 


Bf, = (115) 

where we require the vector to be timelike close to the north and south poles, in the 
physical region — 1 -|- < 0, so that it can be tangent to an observer. If we restrict our 

attention to real helds, and are two spatial vectors in the region — < 0 since 

Efj_v^ = 0 = B^v^. They can however be timelike beyond the velocity of light surface. 
The degeneracy condition on the electromagnetic held is given by 

0 = *(F A F) = = :^E^Bf^, (116) 

where we used F^i, = ^ (2E[^v^] — e^i^^sB'^v^). The second electromagnetic invariant 

/2 = *(F A *F) = = - —{E^E^ - B^B^) (117) 

^ %) ' %) 
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tells us whether the solution is magnetically dominated h < 0 , null /2 = 0 or electrically 
dominated I 2 > 0. The criteria for having a drift velocity of charged particles less than 
the speed of light is to have magnetic dominance I 2 < 0 ; see e.g. 1 ^. 


For the (complex) highest-weight ansatz, the degeneracy condition (116) is written 
as (114). In order to get some intuition it is useful to hrst concentrate on Poincare 


coordinates (f, r, 6 , cj)). The held strength then reads as 


/ 0 {h — 1 ) 01 ( 0 ) —ra[{6) —iqrai{6) \ 


F = 




-h 


0.3(6) 

7r 


-h 


02(6) 


V 


0 o,2i^) — iQ 

0 


03(6) 


(118) 


/ 


It is therefore natural to distinguish four (partially overlapping) qualitative classes of 
solutions type^ 

Poincare magnetic oi = 0 

Poincare electric 02 = 03 = 0 

Poincare nontoroidal ^ a 2 = 0, q = 0 

Poincare generic oi 7 ^ 0 and 02 7 ^ 0 which obey (114). 


A Poincare magnetic solution has no electric held with respect to Any such real 
solution is therefore magnetically dominated. For example, an axisymmetric conhguration 
(g = 0) with real weight h and Oj’s is real and magnetically dominated. (Other real 
solutions can be obtained by superpositions as will be discussed in Sec. 6.1 ) A Poincare 


electric solution has no magnetic held with respect to the 3 + 1 decomposition involving 
the Poincare time t. Any such real solution is therefore electrically dominated. A Poincare 
nontoroidal solution has no components of the electromagnetic held along dcp. This implies 
that the electric held has no toroidal components while the magnetic held (related to the 
dual of F) has no poloidal components. Since the toroidal and poloidal subspaces are 
orthogonal, it is indeed consistent with = 0. In general, there are still toroidal 

magnetic helds and poloidal electric helds but nothing prevents us from canceling one 
such held. The solution can then also be either Poincare magnetic or Poincare electric. 
The Poincare generic solution has no particular electromagnetic property with respect to 
Killing time t. For real helds, there might however be another observer that identihes the 
solution as magnetically or electrically dominated or null. In fact, we will encounter such 
cases in Sec. 


^Poincare magnetic was denoted as Type M in 22 while Poincare electric was denoted as Type E and 


Poincare generic was denoted as Type E-M. Our terminology emphasizes the role of the Poincare time t 
in the 3 + 1 decomposition. 
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Let us mention that this classihcation extends to black hole coordinates. A Poincare 
magnetic solution will also be magnetically dominated in the physical region in the black 
hole patch. Indeed, the vector held 

Vm f 

is real and timelike in the physical region {v ■ v = —sign(l — ^^ 7 ( 0 )^)) and it lies in the 
kernel of F. Observers tangent to that vector held see therefore no electric held. Also, 
a Poincare electric solution admits a real vector in the black hole patch which lies in the 
kernel of icF and which is timelike. The 1-form is given by 

V = (120) 


- 1)/ + /c7(0)V) (119) 


whose norm is —1 everywhere. 


5.6 5'L(2,M) descendants 

Even though the force-free electrodynamics equations are nonlinear, they admit under 
certain circumstances a linear superposition principle. Several sufficient conditions for 
linear superposition were thoroughly discussed in 10,22 . We review some key proposi¬ 


tions here. First, two solutions Ai and A 2 can be linearly superposed if their respective 
currents Ji, J 2 are collinear, A ~ J 2 . Second, a solution A can be linearly superposed 
with its descendant CkA associated with a Killing vector K if the descendant of the cur¬ 
rent is collinear with the current itself, CkJ ~ J- Moreover, under the same assumption 
CkJ ~ J , the gauge held P{Ck)A is also a solution for any polynomial P. While these 
conditions might not be necessary for linearly superposing solutions, they allow us to 
build large classes of solutions. 

In the following, we will check whether solutions with different weight h and U{1) 
charge q can be superposed simply by looking at the current, which will be computed 
in Appendix If they obey the criteria that currents with different values of {h, q) are 
collinear we will be able to linearly superpose them. 

S'L( 2 ,]R) descendants of the highest-weight solution A are dehned from acting with 
Ch_ on A. These gauge helds are not necessarily solutions, except in special cases. From 
the above propositions, descendants of solutions having a current J with components 
along a linear combination of Qq, H_ and /i^ = d6/'y{6) will be solutions because such 
vector helds commute with H_. Another possibility to linearly build a superposition 
of descendants is to start from a highest-weight solution with J ~ ifo- Then since 
[if_, Hq] ~ Hq, Ch_J — J■ We will qualify solutions which obey the property Ch_J cc J 
as “admitting descendants.” 
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5.7 List of solutions 

Here, we list all solutions to force-free electrodynamics with nonvanisliing current that 
we found in our analysis starting from the highest-weight ansatz. We hrst classify the 
solutions according to their highest-weight representation labeled by the (complex) weight 
h and the (integer) t/(l)-charge q and then by their Poincare electromagnetic type. The 
functions Xi{9), i = 1,2, 3,4, 5 obey ODEs in 9 which are described in Appendix [Aj More 
details on the solutions including the held strength, current and Euler potentials can be 
found in Appendix |B| 

(h, q)-eigenstates Two classes of solutions with arbitrary weight h and 17(l)-charge q\ 


Poincare magnetic 

A = 


dh 


■ 2 _ ^g7(l-fcV) 3 


( 121 ) 


where = X^{9-,h,q) and A{h,q) = h{h — 1) -|- k'^q^. The solution is pure 
gauge for h = 0. When g = 0, h 7 ^ 0,1 the solution reduces to (125). 

Poincare generic 




h{h — 1 )X 2 /U^ — kq^X2ii^ + ikq'yX!^^/ 


( 122 ) 


where X 2 = X2{9; A{h, q),Ci = q^). The solution is pure gauge for h = 0. 

(h 7 ^ 0, q = 0)-eigenstates Four classes of axisymmetric solutions with arbitrary weight 
h and one special subcase: 


Poincare generic 


Aikm = cfs'' -A'*-V + A 3 V" ± 


where X 3 = X 3 ( 6 *; h,^), h,^ G C and Ci 7 ^ 0. 

Poincare magnetic 


A(7o) = 


h-1 

TT ,.3 


X4/i^ ± X4 /i' 


where X 4 = X^^iJ)-, A{h)) with A{h) = h(h — 1 ) and C 2 7 ^ 0 . 

Poincare magnetic 

4 = / dht''X,iJ? 


(123) 


(124) 


(125) 


where A, — Xi($; A(h)). 
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• Poincare nontoroidal 


A(,,o) = 


hjji^ ± v^/i^ 


( 126 ) 


where X 2 = X 2 ( 6 ', A(h, g), Ci). 

Poincare electric and nontoroidal - admitting descendants 

A,km = 


( 127 ) 


where X 2 = X 2 { 6 , A(/i, g), 0). It is the special case ci = 0 of (126). 

(h = 0, q 7 ^ 0 )-eigenstates One weight 0 solution with arbitrary 17(1) charge g : 

• Poincare electric 

71(0.9) = (128) 

(h = l,q 7 ^ 0 )-eigenstates One weight 1 solution with arbitrary 17(1) charge g: 

• Poincare electric - admitting descendants 

71 = (129) 


(h = ±ikq, q 7 ^ 0)-eigenstates Two weight ±ikq solutions with arbitrary 17(1) charge 
g: 


Poincare generic 


^{h=±ikq,q) 
^{h=±ikq,q) 




ikqf/ + ig/i^ ± /i^ 


A(h=+ika.n) = <h^A'' k^A + 


( 130 ) 

( 131 ) 


(h = 1 ± ikq, q 7 ^ 0 )-eigenstates One weight l±ikq solution with arbitrary 17(1) charge 
g: 

• Poincare generic - null 

^ A(h{q} = l±ikq,q) = ^ <h^A'^ ^h{q)ai{9)fi^ ± iqai{9)f/ ± 7 a)( 6 ')/i^j. ( 132 ) 

gSZ qSZ 


(h = l,q = 0)-eigenstates Two weight 1 axisymmetric solutions: 
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• Poincare nontoroidal - null 


^( 1 , 0 ) — ® 


± Vc3 + [7ai(6')]V^ 


(133) 


where ai is an arbitrary function. 

• Poincare magnetic 

^(1,0) = ‘h(c2/x^ + Cs/r^). 

(h = 0, q = 0)-eigenstates One weight 0 axisymmetric solution: 

• Poincare electric and nontoroidal - admitting descendants 

^( 0 , 0 ) = 


(134) 


(135) 


This solution is just the S'L(2,]R) invariant solution (62). 


In comparison with 22 , the solutions (124), (126), (130), (131), (132), (133) are new 


and the solutions (128), (129) are given with two branches distinguished by a sign. 


6 Potentially physical near-horizon solutions 

So far, we obtained a list of complex solutions to the force-free equations with various 
types of electromagnetic helds. Physical solutions should obviously be real. In this section, 
we will hrst discuss how to build real solutions from complex solutions. 

Force-free helds are sustained by matter which should travel at subluminal speed or 
at the speed of light in the case of photons. It forces the held strength to be magnetically 
dominated or null. Since the region beyond the velocity of light surface is not physical, 
we will only require that physical solutions be magnetically dominated or null in the po¬ 
lar region between the poles and the velocity of light surface. We will refer to such real 
solutions as the near-horizon solutions which we will list in Sec. 6.2 Finally, solutions 


should have hnite and computable energy and angular momentum hux with respect to 
the asymptotically hat observer and should be solutions to the variational principle. Ac¬ 


cording to the analysis of Sec. 5.3, it leads to the bound (104). The latter set of solutions 


will be refered to as the potentially physical near-horizon solutions, described in Sec. 6.3 


6.1 Reality conditions 

In Poincare and in black hole coordinates the basis /i® is real and therefore real solutions 
are obtained for gauge potentials with real components along the basis functions. Since 
g G Z, h G C we have 




( 136 ) 


27 



























Real solutions can be obtained in various ways. 

First, if h is real and q = 0, then ai{6) is real and the solution is therefore real. 

Second, if the current J and its complex conjugate J* are proportional to each other, 
then one can linearly superpose these solutions to obtain a solution and therefore the real 
and imaginary parts of the gauge field are also solutions to the force-free equations. 

Finally, if the current J and its complex conjugate J* are not proportional to each 
other, one might attempt to find a nonlinear superposition of the solution and its complex 
conjugate with additional nonlinear terms but there is no known systematic procedure to 
do so. 


6.2 List of near-horizon solutions 


Let us now list all real magnetically dominated or null solutions that we could build from 
the complex solutions enumerated in Sec. 5^ At this stage, we list these solutions with 
arbitrary highest-weight h. 


• Nonaxisymmetric, magnetic: 





iq'jjl- , 3 ' 

g2-A(h,g)72 , 


-|- c.c. 


(137) 


where X 5 = h, q) and A(h, q) = h{h — 1) -l- The axisymmetric case g = 0 
is listed below and we have then X^{6, h, 0) = Xi{9, A(h)). 


Nonaxisymmetric, magnetic: 

h{h — l)A 2 /i^ — kq^X2fx^ + ikq'yX^fi^ 
\h*{h* - l)X;fi^ - kq^X;qi^ - ikq-fX;'fi^ 




(138) 


where X 2 = X2{9; A(h, g), ci = g^). The solution is magnetically dominated {I 2 < 0) 
in the range —kq < p < kq. The borderline case /i^ = fc^g^ is a null solution. (This 
is an example of Poincare generic solution which is magnetically dominated.) 


Axisymmetric (magnetic dominance not checked): 


A(h,o)-Ci*h -A3 /i +X^qL ±7/4X3 


(139) 


where h > 1, Ci is real, and ^ > 0. (In the case h = 1, X 3 becomes Xi and the 
solution is not smooth at the velocity of light surface; see Appendix®. It has been 
observed that for h = —1 this solution is magnetically dominated |23| . We did not 
check if it is the case for h > 1. The solution to the nonlinear equation for X3{9) is 
required which we did not obtain here. 












Axisymmetric, magnetic: 


aM _ „h^h 


± X4 /i^ 


(140) 


where X 4 = X4{6]A{h)), C2 is real and arbitrary and h > 2. It is magnetically 
dominated since it is Poincare magnetic. 


Axisymmetric, magnetic: 




dh^^Xi/j^ + c.c. 


(141) 


where Xi = Xi{6, A(h)) and h is complex. For h real we observed that the spectrum 
of h is discrete and the lowest value is greater than 4 in appendix. 


Axisymmetric, magnetic, nontoroidal: 




(142) 


where X2 = X2{6, A{h), ci), h is real. After a numerical check involving X2, it turns 
out that for all Ci > 0 there exists a range of 1 < h < hmax{.Ci) where the solution 
is magnetically dominated for all values of 9. The function hmax tends to 1 in the 
limit Cl —?• 0 and tends to inhnity in the limit ci —)• 00 . It is a solution with no 
toroidal electric field and no poloidal magnetic held. 


Nonaxisymmetric, null: 




hai{6)fi^ ± ±+ c.c. (143) 


EM 

{h=l±.ikq,q) 


Here ai( 6 *) can be complex. We require that ai( 6 *) and a'i{6) vanish at the poles. 
Axisymmetric, null 

Afi% = <F(ai(0)/±7a'i(0)/x3) (144) 

where ai{6) and a'^iO) vanish at the poles but ai{6) is otherwise arbitrary. This is 


a special case of the solution (143) for g = 0. 


6.3 List of potentially physical solutions 

In the following we list the real, magnetically dominated or null force-free solutions with 


nontrivial current and Re(h) = 1 or 1 < Re(h) < 1 and J{6) =0 as dehned in (100), 


which lead to hnite asymptotically hat energy and angular momentum huxes according 


to the analysis in Sec. 5.3 and in 26 . We also do consider linear superpositions. 
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Nonaxisymmetric, magnetic 

We have two classes of solutions from Eq. (137). The hrst one for h{if]) = 1 + i? 7 , 
r] e M: 


A= f 




g2 - ^{h{ri),q)j 




+ c.c. 


(145) 


where X 5 = X^{6] h{ri), q) and A{h{ri), q) = k‘^q^ — rf — When g = 0, the solution 


reduces to (148). 

When Re(h) 7 ^ the constraint J{6) = 0 imposes h = 1. The solution is 


R = ^ <1 

qeZ 


Xs/i' - -X'/i^ 
<1 


+ c.c. 


(146) 


after using (178), where X 2 = X 2 ( 6 ';g^). 


Nonaxisymmetric, null 


The solution (138) has J{6) = 0 for p = ±fcg. This leads to the null solution 

{ikq ± l)X 2 /i^ + igX 2 /i^ + qX^/r^ 


yl(q) = 

where X 2 = X 2 ( 6 '; iifcg, Ci = g^). 

Axisymmetric, magnetic 

A = [ dh^’^Xifx^ + c.c. 


+ c.c. 


(147) 


(148) 


where Xi = Xi( 6 '; A(h)) and h is complex in the range \ < Re(h) < 1. It is not 
clear whether regular solutions exist in that range. Indeed, at least for h real, the 
spectrum of h is discrete for regular solutions and the lowest value is greater than 
4, see Appendix A. Since the solution is axisymmetric, we can check the Znajek’s 
condition 


gee 


(149) 


After taking the near-horizon limit, the generator of the black hole horizon is dt so 
the angular velocity at the Poincare horizon of the near-horizon geometry is VLh = 0. 


We also have = ^{4’) = 0 as shown in (185), and therefore (149) holds. The 
second regularity condition that should be obeyed for extremal black holes only, as 


described in 30 , is also trivially satished. The solution is also regular in the interior 
upon choosing global generators 0 as seen from ([I^. 
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Axisymmetric, magnetic 

The solution (142) has J{6) = 0. For h 

therefore admissible. It reads as 


1, it is magnetically dominated and 


A 




± \/uy 


(150) 


where X2 = ^'2(6', 0, Ci) and Ci > 0. Znajek’s condition does not apply because the 
field is nontoroidal, iq^F = 0. 


• Nonaxisymmetric, null 

Solutions ( |143D with h{q) = 1 ± ikq have J{9) = 0: 


A 


h{q)ai{9)fi^ ± iqai{9)iJ? ± 7%(6')p^ 


+ c.c 


(151) 


We require Oi to vanish at the poles. For g = 0, we obtain the axisymmetric null 
solution: 

• Axisymmetric, null 

A(i,o) = <h(ai(6')/U^ ± 74(6')/!^). (152) 


As discussed in Appendix]^ solutions to the ODEs for Xi, X2 and X5 exist which are 
regular at the north and south poles. However, the functions Xi and X5 are generically 
logarithmically divergent at the velocity of light surface. Since the fate of the velocity of 
light surface is unclear when extending these solutions to the asymptotically flat region, we 
do not exclude them and consider them as potentially physical. A more complete analysis 
of the extension of these solutions to the asymptotically flat region would however be 
necessary to fully settle the issue. 

Generically, highest weight solutions with non-zero weight have a singular field strength 
at the Poincare horizon. However, the energy and angular momentum fluxes are regular 
at the Poincare horizon as a consequence of the restriction (104). One exception is the 
positive branch of (152) which is up to a gauge transformation A(i o) = ai{9)d{t — ^). 
This solution, found originally in 22 , is regular at the future horizon but singular at the 
past horizon. 

The solutions are written in a S'L(2,M) covariant manner and one can choose any 
S'L(2,M) generators related by isomorphisms of the algebra, as discussed around (70). 


7 Summary and conclusion 

This paper provides new analytical solutions to force-free electrodynamics in the near¬ 
horizon region of extremal Kerr black holes. We organized and classihed them according 
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to the highest-weight representation labeled by the complex weight h and the integer 
charge (or angular momentum label) q. Our classihcation rehnes and extends the pre¬ 
vious results in the literature. Since the force-free equations are highly nonlinear, we 
cannot claim to have fully solved them all. Nevertheless, we were able to organize the 
solutions into seven independent and complete classes. Besides such a mathematical clas¬ 
sihcation, we gave a more physical description of the electromagnetic properties of the 
solutions: we scrutinized the linear superpositions leading to real magnetically dominated 
or null solutions which admit hnite energy and angular momentum with respect to the 
asymptotically hat observer. We ended with several families of solutions which, according 
to our mentioned criteria, are the only potentially physical solutions among all the ones 
obtained. However, all of them, except one, are singular at the future Poincare horizon. 

The second main result of this paper is to have developed the consequences of the 
fact that the region beyond the velocity of light surface is not physical, because of the 
lack of a timelike Killing vector held and the consequent impossibility to dehne a vacuum 
for the electron-positron plasma in that region. We discussed that it implies that the 
region near the equator of the near-horizon geometry is unphysical and for this reason 
we only required that physical solutions be magnetically dominated or null in the polar 
regions, namely, between the poles and the velocity of light surface. Our third result is 
more technical and consists in solving numerically the three linear ODEs which appear 
in highest-weight solutions. We showed that they are regular at the poles but pointed to 
logarithmic divergences at the velocity of light surface. 

The main open question left is how to glue the near-horizon solutions to asymptotically 
flat spacetime. We hope to return to this question in the near future. 
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A Relevant ordinary differential equations 

In this section we list the hve ordinary differential equations whose solutions are present 
in Sec. 15.71 

During the resolution of the three coupled nonlinear ODEs in 9 as described in Sec. 
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5.4, we encountered the following nonlinear ODE for X{ 6 ] h, ci, C 2 ): 


X{iX' - 7 X") - 7 (h - l){cl + + X'2) + 727 'XX' (h - l)ci + hkX 


1 2 


+ 7 ' 


(h - l)ci + hkX Uh - 1) hX^ (ci + kX) + [{h - 2)ci + hkX]X'^ 


+X 


{h — l)ci + hkX 


X" \ = 0. 


(153) 


In the case ci 7 ^ 0, one might substitute X = C 1 X 3 and define ^ Then all the 

dependence in ci factors out. We then obtain the differential equation for X 3 ( 0 ;h, 
which is listed below. This ODE was also found in 23 . Upon setting ci = 0 one gets 


another nonlinear ODE. When C 2 7 ^ 0, h 7 ^ 0, C 2 can be factored out of the equation upon 
a rescaling of X. We denote the resulting function as X = C2X\^^{9-^ X{h)). The ODE 
for X4 is listed below. When both ci = C 2 = 0 and h 7 ^ 0 we hnd a linear ODE that we 
denote as X = X\^^{6] X{h)). 

Recall that the unphysical region beyond the velocity of light surface lies in the range 
9^ < 9 < TT — 9^. In the following, we will assume that all functions X{9) together with 
their first derivatives are finite in the physical region 0 < 0 < 0* and n — 9^ < 9 < tt, i.e., 
X{9) < cx), X'{9) < 00 . 

For the extremal Kerr black hole k = 1 and 7, T are given in (|^. In particular it is 
useful to note that ■y{7i — 9) = 7 ( 6 *), 7 ( 0 ) = 7 ( 7 r) = 0, 7 '( 0 ) = 1, 7 '( 7 r) = —1 where 6 * = 0 is 
the north pole and 6 * = vr is the south pole. Also, 9^ = arcsin[\/3 — 1] ~ 0.82 is the lowest 
positive root of 7 ( 6 ^) — 1. All numerical solutions will be plotted only for the extremal 
Kerr black hole. 


1 ) Xi(e;A{h)) 


ODEi[.Yi;A(/i)] = A7 + 


7 ' fc^7^ + 1 


x; + A(h)Xi = 0, 


(154) 


7 / c 2^2 _ ]_ 

where A(h) = h(h — 1). This ODE appears only in solution ( |125 ). From the physical 
requirement that there should not be any singular magnetic flux at the north and south 


poles; see (185): we impose the boundary condition Xi(0) = Xi( 7 r) = 0. 

The equation is invariant under the transformations h —)■ 1 — h which leaves A(h) 
invariant: that is Xi{9-, h) = Xi{9] 1 — h). Moreover, the equation is invariant under the 
reflection 9 ^ n — 9 so that if Xi{9) is a solution, so is Ai( 7 r — 

For the special case A = 0, i.e. h = 0 or h = 1 , the solution is 


Xi(0;O) = C'i + C'2 


5 sin 9 + sin 39 


19 — 16fc2 + 4(3 + Ak^) cos 29 + cos 49 


(155) 


®The function solution to this ODE was denoted as Sh{9) or Sh,m^o{0) in 
solved. 


22 


but was not explicitly 
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However, regularity at the velocity of light surface implies C 2 = 0 and Xi{9, 0) is therefore 
constant which we hx to 0 by the boundary condition Xi(0) = 0. There is therefore no 
solution. 

In general, the differential equations have regular singular points at the zeros of 7 and 
— 1 which are located at Og = 0,9^,% — 0*, tt. Indeed, for each root 9o we have 

' j ' 7^ + 1 1 

-r-= -—— + regular terms. (156) 

7 7 ^ - 1 9 -9o 

Therefore, Frobenius’ method is applicable. In the generic case, the solution reads close to 
the pole 9o as a linear superposition of the power series solutions {9—9o)^^ J 2 n=o ^n{G— 0 oY 
where A± are the two roots of the indicial equation. In case of double roots, a logarithmic 
branch appears. Frobenius’ series converges in the open complex disk that contains only 
one root. 

In the range 0 < 0 < 0* we could start the series from the north pole or the velocity 
of light surface. Now, the indicial roots are 0 and 2 around the north pole while there 
is a double root 0 around the velocity of light surface which leads to a logarithmically 
divergent solution. One might however question whether such a logarithmic divergence is 
admissible. After all, the geometry around the velocity of light surface will be modified 
significantly when considering the asymptotically flat extension of the geometry. If the 
logarithmic divergence at the velocity of light is acceptable, one can simply write an 
expansion close to the north pole Xi ~ 096 *° + 026 ^^ + 0{9'^) with the boundary condition 
Oq = 1 and the choice of normalization 02 = 1. The solution then exists for all (complex) 
values of h and is dehned by a power series expansion. 


Ai(0, A(h)) 


- A « 15.0277 

--- A « 59.1289 
A « 132.4033 
--- A « 234.9028 


0 

Figure 2: Solutions to ODEi[Xi-, A(/i)] for the first discrete values of the parameter 
A(/i) which obey the boundary conditions Ai(0) = t), Xi{6fi) = 1. 



If one removes the logarithmic branch, there is only one constant of integration which 
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we fix by setting Xi{6^,] A{h)) = 1. The solution then takes the form 


X,{9- A{h)) = l + J2 (157) 

n>l 


where the coefficients an{h) can easily be obtained. Now, the physical boundary condition 
Xi(0) = 0 will be obtained only for a discrete spectrum of A{h). The hrst real values of 
A{h) together with the plot of the solution Xi{6-, A{h)) in the physical domain are given 

in Fig. 


2) X2{e-,A{h,q),Ci) 


7 


0DE2[X2; A{h, g), Cl] = X'' + + 


7 


A{h,q) - ^ 
7 ■ 


X2 = 0 


(158) 


where A(h,g) = h(h — l) + k^q^. This equation appears in (147) with ci = q^, in (126) with 


Cl arbitrary and in ( 127[ ) with ci = 0. The equation is invariant under all transformations 
of h, q that leave A(h, q) invariant. In the axisymmetric case g = 0 and when ci 7 ^ 0 we 


require from the analysis around (186) the boundary condition X 2 ( 0 ) = X 2 ( 7 r) = 0 . 


Let us hrst analyze the simplest case ci = 0. This ODE 2 generalizes the one written in 
22 [see their eq. (4.9)] for g arbitrary. We concentrate on solutions which are symmetrical 


around the equator and therefore we only consider the interval 9 = [0, |]. After performing 
the substitution x = sin^( 6 *), ODE 2 takes the form of the generalized Heun’s equation 


„ (« + /? + l)x‘^ - (a + /d + 1 + 0(7 + 5 ) - 5 )x + 07 , ajSx -b v — n 

A-2 / ii/ \ A 2 T / .|\/ ^ 

x{x — l)(x — a) x[x — l)(a: — a) 

where 

a = 2, h = -^, a /3 = -^, a + /? = -^, 7 = 1, ^ = \- ( 159 ) 


Frobenius’ method can be applied. There are poles at the north and south poles and at 
the fake pole x = 2. One could therefore expand in a power series at the north pole and 
it will converge over the region x G [0,1[. It is important to note that the equator x = 1 
( 6 * = |) is not included in the radius of convergence so care should be taken. At the north 
pole we hnd 0 as a double root of the indicial equation. Removing the logarithmically 
divergent branch, we get the following regular convergent power series in the domain 
9e [0,f|: 


A'2(9; A(A,,)) = (160) 

n=0 
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The solution (160) behaves close to the north pole as X 2 ( 0 ) = do where do is arbitrary 


(which we fix to 1 by linearity of the equation) and ^^(O) = 0. The coefficients obey the 
second-order recurrence relation dn+i = Andn + Bndn-i with di = —^do and 


— A 

" ^ 4(n + l)2’ 


Bn = 


2(n- l)(2n-3) - A 
8(n +1)2 ■ 


(161) 


Now, a numerical convergence analysis reveals that the series expansion does not converge 
at the equator a; = 1 unless A = 0. We interpret this by the presence of a source at the 
equator for generic values of A. In Fig. we plot X2{0, A(h, q)) obtained from the series 
expansion truncated to order 20 for some real values of the parameter A{h, q) (the value 
of the function around | for A 7 ^ 0 should be taken with a grain of salt since the series 
does not converge there). 



0 ^ 7T 

2 


0 

10 

20 

40 

80 


Figure 3: Solutions to ODE 2 [X 2 ; A{h,q);t)] for different values of the parameter 
A{h,q). 


Let us now analyze the more general case ci 7 ^ 0. The ODE 2 is still only singular at 
the zeros of 7 , i.e., at points 9 = 0,7i. After performing the substitution 2 ; = sin^( 6 '), the 
ODE 2 takes the following form 


A"(z) + 


z‘^ — Qz + A 
2 z{z-l){z-2) 


A'(z) + 


{z - 2)2ci - AAz 
16z{z — 1 ) 


Xo 


= 0 


(162) 


This makes it clear that z = 0,1 are regular singular points so that Frobenius’s method 
applies. The two solutions of the indicial equation at the north pole are A± = In 

order to avoid oscillations at the poles we enforce ci > 0 from now on. Only the solution 
A+ is admissible since otherwise the solution will diverge at the north pole. In the special 


“In the case Ci = this equation was considered in 22 ; see their (5.23). 


36 














case where A+ — A_ is an integer q which incidentally occnrs for ci = one independent 
solution contains a logarithmic branch which again diverges. Again in this case, only the 
solution which behaves as is admissible. 

In all cases, the regular solution is given for ci > 0 by 


X2{z) = z^+ anz"- 

n=0 


where Oq is an arbitrary constant and 


ai 


02 


A 

— 4A(1 + + 3ci + (2 + Ci)y^ 

32(2 + Cl + 3^) 


For n > 2, we have am+i = Amam + Bmam-i + Cmam -2 where 

^ 6 m{m + i/ci) — A 

4(m + 1) (m + 1 + ycl) ’ 

2A — 4(m — l)(2m — 3) + Ci + 2(9 — 4:m)yA^ 
16(m + l)(m + 1 + ydil) 

(j — _ A _^ 

32(m + l)(m + 1 + i/cl) ’ 


(163) 


(164) 

(165) 


(166) 

(167) 

(168) 


We fix oo = 1 without loss of generality. In the case ci = 0, we recover the recurrence 
relation (160). From (163) we directly see that for all Ci > 0, the function X 2 obeys the 
boundary condition X2(0) = X2(7r) = 0. We again observe numerically that the series 
(163) does not converge at the equator 6* = | unless Ci is fixed as a definite function of 
A(h), 


Cl = Ci(A(h)) (169) 

which asymptotes to 0 for A = 0 and to cx) for A —)• 00 . For example, for ci(lO) ~ 4.90, 
Cl(20) 8.48, Cl(40) 15.02, ci(80) 27.23. In Fig. we plotted the power series 

solution truncated to order 20 to ODE 2 [X 2 ; A(h, g) = 15, ci] for different values of the 
parameter Ci. Note that the boundary condition X2(0) = 0 is true only for ci > 0. 

In Fig. we studied the behavior of solutions (except at 6^ = |) to ODE 2 [X 2 ; A(/i, q), ci 
3] by varying A(h, q) at constant Ci = 3. 
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X 2 (e; A(h)= 15 , Cl) 



3) Xs{e-h,0 


ODEaiXg; K i] ^ - jX'') - j{h - l)(e + ^3^ + X'^) 


+ 7V^3^3 {h-l) + hkX3 


+ r 


{h-l) + hkX3 U/i-1) hX^il + kXs) 


+ [{h - 2) + hkXslX. 


+ X 3 


{h-l) + hkXs 


X'' ^ = 0. 


This ODE appears in solution (123). According to the discussion around (182) we require 
the boundary condition ^ 3 ( 0 ) = X 3 ( 7 r) = 0 . 

When h = 0, the equation reduces to 

de i^deiX^^)) = 0. (170) 

When h = 1, the equation reduces to ODEi[X 3 ( 6 '), A = 0], therefore 

x,{e,h = i,o = x,{e-A = o). (i7i) 

We will not solve this nonlinear ODE here. When h 7^ 0,1, the equation was considered 


in 


22,23 and solved in 23 for the case h =—1. 


4) Xi{9-,A{h)) 

ODE4.Y,iA(ft)]^.Yr+lj3^A'l + A(ft).Y4 + ^-U_.Y," =0 (172) 

where A(/i) = h{h — 1) which we can rewrite as 

A(h) h -2 

ODEi[X4;A(/i)] + ^-y^X4'^ =0. (173) 
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X2(0; A(h), Cl =3) 



A = 0 
A = 10 
A = 20 


A = 40 
A = 80 


Figure 5: Solutions to ODE 2 [X 2 ;A{h,q)]Ci] for ci = 3 and different values of 
A{h,q). 


This ODE appears in solution (124). According to the discussion around (184) we require 
the boundary condition X 4 ( 0 ) = X 4 ( 7 r) = 0 . 

We will not solve this nonlinear ODE here. Note that when h = 2 the equation 
becomes linear with a nonhomogenous term. 


5) X,{e-,h,q) 


" 7 


2h(h- 1)7^ 


(1 — q)'y'^ — q^) 


X' + 


A(h,g)-\ 


T 


X, = 0 


where A{h, q) = h{h — 1) + This linear ODE appears in solution (121). 
It admits the symmetries 


(174) 


A' 5 («; h,q) = X,(e; 1 - h, q) = W(«; K -q) = Xr,{e-, 1 - h, -q). (175) 

We note the special cases 

X5(0;h,g = O) = Xi(0,A(h)), (176) 

X5(0;h = O,g) = X2(0;A = fcV,g'), (177) 

X^{e-h = fq) = X2{9-A = kWq^). (178) 


Since the ODE for Ai was analyzed previously we concentrate on g 7 ^ 0 only. 

There are always two regular singular points in the range 0 < 0 < first ai 6 = 0 
(north pole), and then at 0 = 6 ** = arcsin(\/3 — 1) (velocity of light surface). When 


^This ODE was also found in 
ODE however slightly differs. 


22 


; see their (3.29) where X 5 is denoted as Sh,m- Our analysis of the 
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h is real and for g 7 ^ 0 , there is a regular singular point at the real root of A 7 ^ = 
which is always in the range 0 < 0 < |. There is also an imaginary root which obeys 

sin^ > \J+ \/33) ~ 2.52 (bound reached dX h = q = 1) so it is irrelevant for 
discussing convergence in the interval 0 < sin 6 < 1. 

The two independent solutions behave close to 6^ = 0 as 0'? and 6~'^, while they behave 
close to 6 = 6^ as log(6' — 6**) and {6 — 6**)°. If one only insists in having a solution 
smooth at the north (and south) poles, a solution always exists but it will be generically 
logarithmically divergent at the velocity of light surface. In order to avoid singularities, 
we need to interpolate between the solutions 6*'^ (we assume g > 0) at 0 = 0 and (6* — 0*)° 
at 9 = arcsin(A/3 — 1). This involves a shooting method which will discretize the possible 
values of h as a function of q. We then normalize the solution with = 1. There 

is therefore no more free continuous constant of integration. In the range 0 < h < 1 
and g > 0 the other singularities are not in the range between the north pole and the 
velocity of light surface. For h > 1 and h < 0, they are but we checked that the indicial 
equation around the real pole and imaginary pole has exponents 0 and 2 so an interpolating 
function between the velocity of light surface and the north pole will be smooth. The 
regular solutions for the hrst four real values of h are depicted on Fig. [^for g = 1. 


Xsie, h, q=l) 


- h « 6.0387 

--- h« 9.9507 
h« 13.8111 
--- h« 17.6551 


0 6 ^ 



Figure 6: Solutions to ODE^lX^] h;q = 1] for g = 1 and different values of h such 
that 765 ( 0 ) = 0. 


B Properties of all highest-weight solutions 


In this section, we analyze the properties of all solutions listed in Sec. 5.7 In particular, 
for each vector potential we compute and list its held strength and its current. These 
properties constitute the ID card of each solution. We also compute the canonical Eu¬ 
ler potentials (for dehnitiveness in Poincare coordinates) dehned in Sec. We hnally 
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check for regularity of the solutions at the poles in order to derive the relevant boundary 
conditions for the ODEs that the solutions depend on. 

We recall the definitions A(h) = h{h — 1) and A(h, q) = h{h — 1) + k'^q^. 


{h, g')-eigenstates 

There are two classes of solutions with arbitrary nonvanishing highest-weight h and 17(1)- 
charge q, describing stationary and nonaxisymmetric field configurations. 


Poincare magnetic 




2 m(l - ^V) 




q^ — A(h, 5 ) 7 ^ 


— A(h, q)^"^ 


^ — iq{l — k‘^'y^)w^ + (1 ~ h)k'y‘^w‘^ — 


+ X5{w^-w^) 


J(7,) = -^^X'^hk 


X.+ 


2(1 - h)7y 


(1 — k‘^'y‘^)(q‘^ — A(h, g) 7 ^) 




$77+ 

yfp 


where A 5 = X^{9]h,q). Since the current keeps its direction upon changing h, g, one 
might linearly superpose solutions with different h, q. The current J{h,q) and its complex 
conjugate {J(h,q))* are proportional to each other. Then, we can obtain a real solution by 
adding up the vector potential A(^h,q) and its complex conjugate. 

We do not consider the trivial case h = 0, because it gives no field strength F = 0. For 
h = 1 , g y 0 the solution exists and the function X 5 becomes X^{9; 1 , g) = X2{9, k'^q^, g^). 


In the case g = 0, h y 0,1, we recover the Poincare magnetic solution (125) after using 
the identity X^{9-, h, 0 ) = Xi{9; h). 

In order to get physical insight, it is useful to derive the functional expression of Euler 
potentials. We only consider Poincare coordinates and g 7 ^ 0. The held strength describes 
a stationary and Qo-eigenstate conhguration with = 0 and df = 0 according to the 
analysis of Sec. Therefore it takes the form (57) with 


V’(a^) = 


de'ip2{r,9) = -iq\ 


1 - ^272 




g2-A(h,g)72 /X 5 


(179) 


The Euler potential y is singular at the Poincare horizon r = 0, unless h is nega¬ 
tive. In the special case h = 1 and g y 0, we can easily integrate ^2 and obtain 
'^2{9) = ln(A 5 ) const. 
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Poincare generic 


A 


{h,q) 


= 


h{h — l)X 2 /i^ — /cg^X 2 /i^ + ikq'yX2fi^ 




X 2 


'yX2{ikqw^ — {h — l)t(;^) + 


k\^ + {h-lY 4 


k 


+ kq'^w^ + i{h — l)gw® j 


hk 


Jih,q) = $"A«—X 2 


{h - lyYQo + kq- 


y/2T 


iq{h — 1 ) 1 / 


where X 2 = X2{9-, A(h, q),ci = q^). 

The case h = 0 is trivial since F = 0. The current is not proportional to its complex 
conjugate, unless for the special cases g = 0orh = lorh = l + z/i for any real /i. The 


case g = 0 coincides with the solution (127) and it will be analyzed below. The case h = 1 


actually coincides with the solution (121) for h = 1. It was just analyzed in the previous 


subsection. The third class is an independent real solution. 

In the generic case g 7 ^ 0, h 7 ^ 0,1, the held strength describes a stationary and 
Qo-eigenstate conhguration with ig^F 7 ^ 0. The canonical Euler potentials in Poincare 


coordinates can therefore be written as (46) where 


Xi{r,e) = h{h - 1 ) 


X 2 

r>h—l ' 


X2{r) = 


ikq 1 


= 0 


h — Ir' 


(180) 


[h 0, q = 0)-eigenstates 

Poincare generic 


F(,,o) = c?4>"-'X3"-2 


Cl 7 ^ 0 


h 


hXs + 


Fh^/^Xsw^ + 7((h - 1) + hkXyX'^w^ - —X'^X^w^F 

kj 

(h — 1 ) \ ^ ^4 ^ v2„„5 


X^w^ - hX^w^ 






{h-l)A{e;h,OQo-hXs 


h^/T-A{e;h,0 $ 77 + 

7^ y2r 


h{h-l)^^ 3 ^ h(h - 1) V? ,,4 




7r 
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where A{6] h, is given by 


{T[y\h - 1 + hkX^f - h^Xl])A{e; h, 0 = 

- 1 + hkX3)][-hX3 + k-f\h - 1 + hkXs)] 

+ 2h'yXs{h - 1 + hkX^)iX'^ - h{h - l)^‘^X'i (181) 


where X 3 = X^{6]h,^). In the case h, ci G M and ^ > 0, the expression of the vector 


potential is real. For h = 1, we get the Poincare generic solution (134) [because 

X 3 ( 6 ', 1 ,^) =Xi( 6 *; 0 ) = constant]. 

The Euler potentials for this stationary and axisymmetric conhguration fall in the 


category (28) where in Poincare coordinates 


'ip{r,e) = 



= fl (?/>) = Cl 


h-l 


h 




(182) 


Since 2TT'i(j{r, 6) is the magnetic flux through the loop of revolution defined by (r, 6), the 
requirement of having no singular magnetic flux at the north and south poles is equivalent 
to the boundary conditions X 3 ( 0 ) = X 3 ( 7 r) = 0 . 


These potentials allow us to recognize the solution as the one described in 23 upon 
identifying their quantities in terms of ours as 


a = -h, m^{ciX3)\ = = ° = 

Poincare magnetic 


Aihn^ — Co<F^ 


h/i,o) — >-2' 

— >-2' 




F(hm = ^ ±hX^ ^ + k'jX^w^ — - — X'^w^ + hX^iw^ — w^) 


1 

kj 


C(9-h) _ 

Jih,o) = 4^ h{h - 1 )X 4 Qo - kX^ .xJj. .1 T h{h - 1 )X 4 p 


-1 

TT ,.3 


7(A;V - 1) V2T 


7 


where C{6] h) is given by 


C{e- h) = hy [{h - l)Xf/^ - 1 ) - 


. 2^,3 


(183) 
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Since dtF = 0, the Euler potentials for this stationary and axisymmetric conhguration 
are given by (28) where in Poincare coordinates 





/(V’) = = 0 


(184) 


Since 6) is the magnetic flux through the loop of revolution dehned by (r, 0), the 

requirement of having no singular magnetic flux at the north and south poles is equivalent 
to the boundary conditions X4(0) = X4(7r) = 0. In turn, regularity of the current at the 
poles then requires h> 2. 


Poincare magnetic 


0 ) = 


F{h,o) — 




h-l 


k'yX.nP' — - — X,vj^ + hXiiw'^ — 
k'j 

^_^2kiX[ + hk'yXiil - 


w 


-1 + A;V)7r 


V2T 


where Xi = Xi(6'; A(h)). Currents with different values of h are collinear so one might 
linearly superpose such solutions. 

The Euler potentials for this stationary and axisymmetric configuration are 




(185) 


Since 27r'0(r, 6) is the magnetic flux through the loop of revolution dehned by (r, 0), the 
requirement of having no singular magnetic hux at the north and south poles is equivalent 
to the boundary conditions Xi(0) = Xi(7r) = 0. 

Poincare nontoroidal 




\h-l V I 


F(/j,0) = ^Xs I - l—w^ - 


^ 2.,.2 (^^^4 
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hh^ 

^2r 


<i>Er4 
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^/2V 
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where X 2 = X 2 { 6 -, h, ci). 

For h = 0 this solution has vanishing held strength so is pure gauge. For ci = 0 we 
get a Poincare electric and nontoroidal solution, which has the special property to admit 
descendants solutions. The solution is real in Poincare coordinates for ci > 0, h G M. 

Since i^F = 0, the electromagnetic held in terms of Euler potentials takes the special 
form (32) where in Poincare coordinates, 


X{r,9) = h 


X 2 

yh —1 ' 


Hx) = =F\/nx 


(186) 


We observe that in order to prevent singular line currents we need to enforce a vanishing 
polar current I at north and south poles, which requires the existence of the boundary 
conditions df 2 ( 0 ) = X 2 { 7 i) = 0 for ci 7 ^ 0. For ci = 0 the polar current vanishes but the 
electrostatic potential y is constant on the north and south poles. 


{h = 0, q ^ 0 )-eigenstates 

Poincare electric 


A \q f —dO 1 

^(0,9) = A^e ^ /X , 

F(o,q) = + iqw^^ 


Ji 


(O.q) 


±f3-d9 
= A^—e •' -y 


-{q^ + 7^)Qo + fg(±g/i^ - /i^) 


The held strength is singular at either the north or south pole depending upon the sign. 
The solution might however be interesting if it is split at the equator with regular north 
and south branches. 


In terms of Euler potentials and in Poincare coordinates, we are in the case (46) where 

(187) 


yi(r,6') = re 


±( 3-de 

- T' J q' 


72 ( 76 ') = 0 , k( 0 i) = 0 


The vanishing of X 2 is related to the absence of magnetic held. 


(/i = 1 , g 76 0 )-eigenstates 

Poincare electric - admitting descendants 

A rTk \ <7 i r ~ 1 
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The field strength is again singular at either the north or south pole depending upon 
the sign. The solution might however be interesting if it is split at the equator with 
regular north and south branches. The direction of the current does not depend upon q 
and therefore one can linearly superpose solutions with different g’s. The solution admits 
descendants. 


In terms of Euler potentials and in Poincare coordinates, we are again in the case (46) 
where 


Xi{0) = e 


±f^de 

- J 7 


X 2 {r, 0 ) = 0 , K{(j)i) = 0 


(188) 


{h{q) = zhz/cg, q ^ 0)-eigenstates 

The two following classes of solutions feature a charge-dependent weight h. 


Poincare generic 

A{h{q),q) = ^ + iq^x^ S 2 /i^ j, S 2 = -1 or 1 

F(h(q),q) = ^ ^±S2W^ — S 2 w‘^ + =F iqw^ + iqw^^ 


JiKqU) = 




[Tq + - 1 ) - Ki ± kq)Y]Qo 




7 


The solution is pure gauge when g = 0. In terms of Euler potentials and in Poincare 


coordinates, we are in the case (46) where 


Xi(p^) = T , y2(r) = T-, fi:(0i) = O 


(189) 


Poincare generic 


^(h(q),q) — kfi^ + 


F[h[q),q) = + ikq{Tw^ + w®)) 

Jwq),q) = $"A^^(±g + zfcf ) ( (z ± kq)Qo + 


q (^^H+ 4 

7^ V y2r ^ 
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In terms of Euler potentials and in Poincare coordinates, we are in the case (46) where 


, . k , , 1 


X 2 (r)-±^, 

k(0i) = 0 


(190) 


{h{q) = 1 zh ikq, q ^ 0)-eigenstates 

Poincare generic - nnll 


A 






hai{6)iA ± iqai{6)jj^ ± 'ya\{6)jj? 


\h(q),q) 




qa\{ 6 ){±.w^ — 'uA) + iqai{ 6 ){'^w^ + w^) 
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{Hq),q) 




(1 ± ikq)[q^ai - 'ydei'ya'^)] 

yr 


/c7^Qo - 


^H+ 

v/2r 




where a\{9) is an arbitrary function. It is a null solution = 0). The current is 

/ 2 

nonvanishing for g 7 ^ 0. Indeed, the current vanishes when ai{ 6 ) obeys a'/+—^ai = 0. 
After a closer look at this differential equation, we conclude that a solution is given by 
X 2 { 9 ; A{h{q),q) = 0,Ci = g^). The constraint A(h(g),g) = 0 however implies g = 0 in 
contradiction to our assumption g 7 ^ 0. 

In terms of Euler potentials and in Poincare coordinates, we are in the case (46) where 

ha{9) 


Xi{r,9) = 


r^h—1 


X2{r) = ±-, k ( 0 i ) = 0 

r 


(191) 


[h = 1, q = 0)-eigenstates 

Poincare nontoroidal - null 


ai{9)ij,^ ± a/c 3 + [ya'i 




^(1,0) — “f* 

F(1,0) = ±\/cs + [701(6')]2^;^ - 'ya[{ 9 )w'^ 




deixa'i) 


-kqQo + 


•hi/. 


± 


a'l 


V2r y/c 3 + [ 7 a) 


r/^ 


It is a null solution {F^^F^^'' = 0). This class of solutions does not overlap with the class 
above. The current is vanishing when d 0 {'^a'^ = 0, i.e., when ai = C 2 + / ^d9, where Ci 
and C 2 are two real constants. Regularity Exes Ci = 0 so only constant ai solutions obey 
Maxwell’s equations. 

In terms of Euler potentials, the held strength reads as (32) where 


X(r,0) = ai(0), /(x) = T^cs + [qx' 


(192) 


47 































We see that we need C 3 = 0 in order to have a regular conhguration (no polar current on 
the 9 = 0 axis). We also impose that ai{9) must be regular at the poles. 


{h = 0,q = 0)-eigenstates 

Poincare generic - admitting descendants 



The solution with C 3 7 ^ 0 is singular at the poles. Indeed, 

f dO cos( 6 *) 

J 2 


In 


, 9 

tan I - 


(193) 


Therefore we fix C 3 = 0. The solution then becomes electric without toroidal fields. 
In fact, it is just the maximally symmetric solution. It is related to (135) by a gauge 
transformation. 
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